Single transverse spin asymmetry of forward neutrons 
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We calculate the single transverse spin asymmetry Aiv(t), for inclusive neutron production in 
pp collisions at forward rapidities relative to the polarized proton in the energy range of RHIC. 
Absorptive corrections to the pion pole generate a relative phase between the spin-flip and non-flip 
amplitudes, leading to a transverse spin asymmetry which is found to be far too small to explain 
the magnitude of An observed in the PHENIX experiment. A larger contribution, which does 
not vanish at high energies, comes from the interference of pion and ai-Reggeon exchanges. The 
unnatural parity of oi guarantees a substantial phase shift, although the magnitude is strongly 
suppressed by the smallness of diffractive Tip — > aip cross section. We replace the Regge ai pole by 
the Regge cut corresponding to the np exchange in the 1^ S state. The production of such a state, 
which we treat as an effective pole a, forms a narrow peak in the Stt invariant mass distribution in 
diffractive irp interactions . The cross section is large, so one can assume that this state saturates 
the spectral function of the axial current and we can determine its coupling to nucleons via the 
PCAC Goldberger-Treiman relation and the second Weinberg sum rule. The numerical results of 
the parameter- free calculation of Am are in excellent agreement with the PHENIX data. 

PACS numbers: 13.85. Ni, ll.80.Cr, 11.80.Gw, 13.88.-|-e 



I. INTRODUCTION 

The single transverse spin asymmetry of neutrons was 
measured recently by the PHENIX experiment at RHIC 
[l| in pp collisions at energies y/s ~ 62, 200 and 500 GeV. 
The measurements were performed with a transversely 
polarized proton beam and the neutron was detected at 
very forward and backward rapidities relative to the po- 
larized beam. Preliminary results are depicted in Fig. [1] 
An appreciable single transverse spin asymmetry was 
found in events with large fractional neutron momenta z. 
The data agree with a linear dependence on the neutron 
transverse momentum qx, and different energy match 
well, what indicates at an energy independent Aff{qT). 

Usually polarization data are more sensitive to the 
mechanisms of reactions than the cross section. Below we 
demonstrate that the large magnitude of the single trans- 
verse spin asymmetry of forward neutrons discovered in 
reveals a new important mechanism of neutron pro- 
duction ignored in all previous studies of the reaction 
cross section. 

At the same time, neutrons produced with xp < 
show a small asymmetry, consistent with zero. This fact 
is explained by the so called Abarbanel- Gross theorem 
which predicts zero transverse-spin asymmetry for parti- 
cles produced in the fragmentation region of an unpolar- 
ized beam. This statement was proven within the Regge 
pole model illustrated in Fig. [31 The amplitude of the re- 
action pt+p X + n squared, Fig. [2^, is related by the 
optical theorem with the triple-Regge graph in Fig. [^Jj. 
According to Regge factorization the proton spin can cor- 
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FIG. 1: (Color online) Single transverse spin asymmetry An 
in the reaction pp — > nX, measured at — 62, 200, 500 GeV 
[Jl (preliminary data). The asterisks show the result of our 
calculation, Eq. (|38|l . which was done point by point, since 
each experimental point has a specific value of z (see Table[l|. 



relates only with the vector product, [kx k'] , of the proton 
momenta in the two conjugated amplitudes, as is shown 
in Fig. ^p. According to the optical theorem these mo- 
menta are equal, fc = fc', so no transverse spin correlation 
is possible. Regge cuts shown in Fig. [2]: breakdown this 
statement, but the magnitude of the gained single-spin 
asymmetry calculated in [3] , turned out to be very small, 
less than 1%. 
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FIG. 2: Graphical representation of the cross section of in- 
clusive reaction pt+P X + n with a polarized proton beam 
and neutron produced in the fragmentation region of the un- 
polarized proton, and the related triple Regge graph (a-b). 
The Regge cut correction to the previous graph (c). 



II. SPIN STRUCTURE OF THE PION POLE 

The importance of the pion pole for neutron produc- 
tion has been known since the pioneering paper [4]. Its 
Reggeized version for leading neutrons in DIS was pro- 
posed in A further development, which includes 
the spin structure of the amplitudes was published in 
Q- Contrary to the conventional wisdom that the pion- 
nucleon vertex is pure spin-flip, a large non-flip term in 
the amplitude appears if z < 1. In Born approxima- 
tion the pion exchange in neutron production, depicted 
in Fig. [3l in the leading order in the small parameter 
rriN / \/s has the form [3], 



(T3 qL + (J ■ qr 
\ z 



i^<f{qT.z), (1) 



where u are the Pauli matrices; fp^„ are the proton or 
neutron spinors; qr and 



1 - z 



IL 



(2) 



are the transverse and longitudinal components of the 
momentum transfer respectively. 
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FIG. 3: Graphical representation of the cross section of in- 
clusive neutron production in hadron-proton collisions, in the 
fragmentation region of the proton. 



At large z thepseudoscalar amplitude (j>^{qT,z) has 
the Regge form 



X A,+p^x(A4), (3) 



where = (1 — z)s; the 4-momentum transfer squared 
(t) has the form. 



and ?/7r(0 is the phase (signature) factor, 

The second term contains the pion pole, 

2 1 
Rer]^{t) « — - - 



t 



(4) 
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In what follows we assume that the pion Regge trajectory 
is a linear function of t, aTr{t) = a'^(t — m^), where a'^ « 
0.9 GeV-^ 

The effective vertex function is parametrized as, 



Rit 



where the pion-nucleon coupling 5?+„„/87r = 13.85. The 

t-slope parameter i?J incorporates the t-dependences of 
the coupling and of the irN inelastic amplitude. Al- 
though it is not well known, its value is not really impor- 
tant for us, since we concentrate on the small t region. 
For further calculations we fix i?^ — 4GeV^^, which is 
naturally related to the nucleon size. 

The amplitude H]) is normalized as Mj^al^i^ = 
X] l^7r+p-*jf Correspondingly, the differential 

X 

cross section of inclusive neutron production reads. 



dz dq^ 



= \t\a 



(8) 



This pure pion pole model has two obvious shortcom- 
ings: (i) the cross section Eq. (|8]) substantially overshoots 
data [7[; (ii) no single transverse spin asymmetry is pos- 
sible because the spin-flip and non-flip terms in the am- 
plitude (dJ have no phase shift. 



III. SPIN EFFECTS GENERATED BY 
INITIAL/FINAL STATE INTERACTIONS 

It is known that initial and final state interactions 
(ISI/FSI) can generate spin effects [S], if they are absent 
in the Born approximation. In the process under discus- 
sion ISI/FSI essentially modify the spin amplitudes of the 
process, suppressing the magnitude of the cross section 
and generating a nonzero relative phase shift. 

Absorptive corrections look simpler and factorize in 
impact parameter. Therefore, we perform a Fourier 
transform on the amplitude ([1]) to impact parameter rep- 
resentation, multiply it by the survival probability am- 
plitude S{b, z). 



fp^nib, z) = fp^nib, z) S{b, z). 



(9) 
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and then return back to the momentum representation where 
The Born amphtude ([IJ in impact parameter reads, 



<j3qLei{b,z)-i^ef{b,z) 



(10) 



J 



'^{b,z) ^ J d'qe^'^^^^^iqT,z) = N^iz)lt^Ko{b/(3) + ^—^[Ko{eb)-K„{b/(3)] 



'^^{b,z) = - I d?qe^'^^{b■q)cp^{qT.z)^Nl{z)l^^K,{b/P) + ^-^ 



eK,{eb)--K,{b/p) 
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and 



X) 



z{ql + ml), 

\ [Rl - «; ln(l - z)\ 
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As is illustrated in Fig. [31 in the rapidity interval cov- 
ered by the pion, Ay — \ ln(l — z)|, no particles are pro- 
duced. Since inelastic ISI/FSI of the spectator partons 
tend to fill this rapidity gap, its survival probability leads 
to a suppression. It was found in Q that the partonic 
system, which participates in ISI/FSI with the target nu- 
cleon, is a 5-quark color octet-octet dipole. The ampli- 
tude suppression factor S'(6, z) was found in Q to be 
rather small, leading to a strong reduction of the cross 
section of leading neutron production. The magnitude of 
the partial cross section, which is expected to be univer- 
sal, was found to agree well with ZEUS data for neutron 
produced in DIS 8], and with NA49 data for colli- 
sions 3 • Although these data confirm the shape of the z 
distribution measured earlier at ISR [l3| , the latter over- 
estimated the normalization. 

Following ([9]) wc applied the absorption factor to the 
Born spin amplitudes Eq. ()10p . and performing Fourier 
transfer back to momentum representation, arrive at. 
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linUQT,z) = 4^1^ / dbbJiibqr) 



X Ki{b/f3)S{b,z). 



(19) 



We are now in a position to calculate the single trans- 
verse spin asymmetry. 



Ajv((7T, z) 



where 



2qTqLVzJ2x l'/'o(gT,^)| \(t)s{qT,z)\ 
^<1lJ2x \<l>o{(lT,z)f + qTJ2x l0s(9T,^)|^ 



X sin((5s - So) , 



tan 5q , 



Im (f>o,s{qT, z) 
Rc (l)o^siqT,z) 



(20) 



(21) 



The results are depicted in Fig.|4l It is clear that our cal- 
culations significantly underestimate the PHENIX data 
in the range of qx covered by the experiment. 
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FIG. 4: (Color online) Single transverse spin asymmetry of 
leading neutrons related to single pion exchange corrected for 
absorptive corrections, as function of qr- The curves from 
bottom to top correspond to z — 0.5, 0.7 and 0.9. 



IV. OTHER SOURCES OF SPIN EFFECTS: 
ai REGGEON? 

In addition to pion exchange, other Regge poles R — 
p, a2, ai, etc. and Regge cuts can contribute to 
the pp — 7> nX reaction. We should remind that the total 
collision energy is shared by the reaction tt + p — >■ i? + p 
depicted in the upper blob of the diagram Fig. ([S]), and by 
the rapidity gap, rather unequally: while the former can 
be very large, = (1 — z)s at high energies, the cm. 
energy for pion exchange is rather small, s' = sq/{\ — z), 
and independent of s. 
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FIG. 5: Graphical representation for the interference between 
the amplitudes with pion and Reggeon exchanges. 



at high energies. 

The aiNN vertex is known to be pure non spin-flip 
[ill [l^ . Therefore, it should be added to the first term 
in Eq. ([T5|). 



2mNqL <at c 

(22) 

where the longitudinal polarization vector of ai reads 

[a, 



t, 0, 0, go 



and the transferred energy 



qq = Ep- En = qL + 0{mN/Vs) 



(23) 



(24) 



In the Born approximation 
a' 



it) ria^it) 



X (l-^)-"-^'^^a>^x(^^. 



and 



tg 



(25) 
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The amplitude (j22p contains three unknowns, to be 
fixed before numerical evaluation: 



• The amplitude ^a+p^xi-^'^. 

• The oi-nucleon vertex G + 

J- ajp 

• The Regge trajectory (t) 



i); 



Notice that the general structure of the amplitude 
Eq. (pSj) is valid for any axial-vector state. In what fol- 
lows we find the ai pole to be quite a weak singularity, 
and conclude that the spectral function of the axial cur- 
rent is dominated by the contribution of np pair in the 
1~^S state. So we will replace the weak ai pole by an 
effective singularity a. 



Summing over different produced states X and using 
completeness one arrives at the imaginary part of the am- 
plitude of the process n+p^R + p at cm. energy M^. 
The production of natural parity states, like p, 02, etc. 
can proceed only via reggeon exchange, therefore these 
amplitudes are strongly suppressed at RHIC energies by 
a power of Mx (dependent on the Regge intercept) and 
can be safely neglected everywhere, except at the region 
of very small (1 — z) ^ so/s, unreachable experimentally. 

Only the unnatural parity states, which can be diffrac- 
tively produced by a pion, like the ai meson, or p-n in 
the axial vector or pseudo-scalar states, contribute to the 
interference term in the neutron production cross section 



A. Diffractive production of ai meson and 
non-resonant np{l~^S) pairs 

The amplitude A^+^^^{M'^) is normalized as, 

X ^ 

X \/ d(T{Tip aip)/dpTlpr=o (27) 

Evaluation of the ratio of the forward diffractive-to- 
elastic cross sections in this expression is a subtle prob- 
lem. The cross section of diffractive ai production is so 
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small that it has been escaping detection for a long time. 
Eventually ai production was observed in n+p — Stt+p 
performing a phase shift analysis of high-statistics data 
[1^ [T5|. The ai resonance was only detected in this pro- 
cess as a fast variation of the relative phase of l+S*. 

The p7r(l+S') mass distribution depicted in Fig. [S] 
forms a strong and narrow peak, related mainly to the 
Deck mechanism fl6|. The data combined for lab 
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FIG. 6: (Color online) Variation of I'^S wave intensity with 
Stt invariant mass, for < 0.05 GeV'^. Data for Eiab ~ 
63 and 94 GeV are summed The curve is fitted with 

parametrization (1281) . 

energies 63 and 94 GeV were taken with the small-gr 
trigger. We normalize this distribution on the mea- 
sured cross section in the invariant mass interval 1.0 < 
< 1.2 GeV which is 54 ± 8 ^ib and 48 ± 7/i6 at 63 
and 94 GeV respectively. With the measured p^-slope 
B = 14.3 ± 0.17 GeV~^ we can evaluate the forward 
cross section for this mass interval at da/dp^]^^ ^ — 

0.73 ± 0.07mb/GeV^ averaged for the two energies. 

In order to normalize the mass distribution plotted in 
Fig. [6] to this cross section, we fitted the data with a 
simple Gaussian parametrization. 



d<j 



dM^T^ dq^ 



qT=0 



= X exp 



(28) 



with parameters = (1.117±0.004) GeV, F^ = (0.226± 
0.004) GeV, and A = 3.91 ± 0.03mb/GeV^ We also 
made a 10% correction for the nonzero values of < 
0.05 GeV^ of data in Fig. [HI Notice that the central mass 
of the peak is significantly smaller than TOai = 1.28 GeV 

Eventually, integrating the normalized distribution 
over the invariant 37r mass under the a-peak in Fig. [5] 
we arrive at the 7rjo(l+5') forward production cross sec- 
tion. 



da^irp — )■ npp) 



dp^ 



1 + S,pr=0 



1.67mb/GeV2 



(29) 



Now we have to extrapolate the diffractive cross sec- 
tion Eq. ((29|) , from the energy range of the fixed target 



experiments [15|, si ^ 119— 177 GeV'^ up to the substan- 
tially higher cm. energies M'jr = (1 — z)s of RHIC. We 
extrapolate the single diffractive cross section with the 
following energy dependence. 



d(Jsd{M'j^) _ dasdisi) 



dq^ 



dq'^ 



Ml 

Sl 



0.16 



(30) 



where K{s) is the survival probability of a large rapidity 
gap , which in the eikonal approximation has the form 



K^p{s) = 1- 



1 



^ BYAs)+2B:1{s) 



7rp 



+ 



(47r)2 SJ(.) [BZ{s) 



5J| 



(31) 



For the energy dependent total irp cross section a^^^ [s) 
we rely on the parametrization given in [l7| . For the 
elastic and single diffractive slopes we use B'^p{s) = 
Bq + 2a']p ln(s/so), where Bq ~ 6GcV^^ for elastic and 
Bq = 9 GeV~^ for single diffractive processes, a'p = 
0.25 GeV"^ 

Notice that the replacement of the oi pole in the dis- 
persion relation by the effective singularity a also helps 
to settle several problems related to the assumed ai dom- 
inance for the axial current. In particular, this assump- 
tion leads to so called Piketty-Stodolsky paradox [l3l.[l9|. 
Namely, based on the PCAC hypothesis, which relates 
the pion pole and heavier states contributions in the dis- 
persion relation for the axial current, one arrives at the 
conclusion that diffractive tt ai and elastic tt — > tt 
cross sections should be equal, while experimentally they 
are different by more than one order of magnitude. Re- 
placing the oi pole by the effective singularity a solves 
the problem ^[lO-Sl. 

Summarizing, we found that the production cross sec- 
tion for the ai meson is too small to produce a sizable 
contribution to neutron production. It should be replace 
by the more significant production of a irp I'^S state, 
which forms a narrow resonance-like peak in the 37r in- 
variant mass distribution. Thus, we introduce an effec- 
tive "pole" a in the dispersion relation for the axial- vector 
current, and predict its production cross section in np 
collisions at high energies. 



B. aNN coupling 



Like for the pion we parametrize the oi-nucleon vertex 
in Eq. ^ as Ga+pn{t) = 9a+pnexp{Rlt). The slope 
parameter is poorly known, however in the small t 
region under consideration it is not of great importance. 
Like for the pion vertex, it is natural to expect the slope 
to be related to the nucleon size, so we fix — R^ — 
4 GeV^ for further calculations. 
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The a-nucleon coupUng Qa+pn can be estimated based 
on PCAC. Although it is tempting to interpret the 
Goldberger-Treiman relation and Adler theorem as pion 
pole dominance, the pion pole does not contribute in 
either /3-decay or high-energy neutrino interactions, be- 
cause of conservation of the lepton current (neglecting 
the lepton mass) [H, El, Ull. In order to have PCAC 
heavy axial states contributing to the dispersion relation 
for the amplitude of the process must miraculously repro- 
duce the pion pole. If we replace the combined contribu- 
tion of the heavy state by an affective pole a [H, f^lj , 
the Goldberger-Treiman relation relates the a and pion 
poles. 



"\/2/a gaNN /tt 9ttNN 



WIN 



(32) 



In the second Weinberg sum rule the spectral func- 
tion of the vector current can be represented by the p- 
meson pole. Correspondingly, the axial spectral function 
is saturated by the a meson, because the pion does not 
contribute to the second Weinberg sum rule. Then one 
arrives at the relation. 



/a /p 



\/2m 



7p 



(33) 



where 7p is the universal coupling (pNN, ptttt, etc). 



= 2.4. 



Thus, for the a to pion couplings ratio we get. 



girNN 2mjM fp 



0.5. 



(34) 



A different expression was derived in [1^ , based on the 
first Weinberg sum rule, which includes the pion pole. 



gaNN 



mampf^ 



g^NN 2mN{fp - m2/2)i/2 ' 



(35) 



Numerically the results (|34|) and (l35t are very close in 
the chiral limit of massless pion. Otherwise they differ 
by 15%. 



C. Regge trajectory of the "a-pole" 

So far the narrow a-peak in the spectral function of 
the axial current could be treated as an effective pole 
replacing the real one oi, which was found to be too 
weak. However, the Regge singularity in the complex 
angular momentum plane, related to the ir-p exchange, 
is a Regge cut rather than a pole. The trajectory of the 
cut can be expressed in terms of the tt and p Reggeons, 



i^p{t) ^ a^{0) + ap{0) - 1 + 



a'^a'p 



t. 



(36) 



For further numerical evaluations we fix q;(0) = 1/2 and 
so Q;^p(i) = -0.5-f 0.45t. 



a'p = 0.9 GeV" 



Correspondingly, the phase (signature) factor for the 
unnatural parity a-Reggeon exchange reads, 



ria{t) = -i-tg [Kaa{t)/2] 



(37) 



where aa{t) — a.,rp[t). 

This factor provides a significant phase shift = 7r/4 
relative to the pion, and this phase shifts rises with t 
up to the maximal value of7r/2att = — 1 GeV^ . This 
interference looks like a promising source of a significant 
single transverse spin asymmetry. 

Notice that the intercept of the Tr-p cut turns out to 
be rather close to the intercept of the ai Regge pole, 
aai(O) = —0.43, which corresponds to a straight Regge 
trajectory with the universal slope crossing the position 
of the oi pole on the Chew-Frautschi plot. 



V. SINGLE TRANSVERSE SPIN ASYMMETRY 
GENERATED BY -n-a INTERFERENCE 

Eventually, we are in a position to perform a parameter 
free calculation of the o-tt interference contribution to 
single transverse spin asymmetry of neutron production, 

AN\qT,z) — qT . . , (1 - z) ; —2 



ffg+pn / dcr.^p^ap{Mjr) / dp^ 

57r+pn \dcr^p^T,p{M]^)/dp'^\ 

Pt=0 



1/2 



(38) 



where Aa{t) — a^(t) — aa{t). 

The specific value of z is important for the factors that 
appear in ([38]), in particular for t given by (l2])-(|4]). We 
extracted values of z for each experimental point shown 
in Fig.[T]) from data [1| as is described in |Appendix A| and 
presented in Table H] The results of the calculation with 
Eq. (|38l) are plotted in Fig. [1] by asterisks. The agreement 
with data is excellent, but we estimate the uncertainty 
of the normalization of the theoretical results at about 
30%. 

Notice that in these calculations we neglected the ab- 
sorptive corrections, because they have very little influ- 
ence on An, Eq. ([55]) . Indeed, as we have seen above, 
they cause a tiny variation of the relative phase, which 
can be neglected. It was also demonstrated in 7] that the 
spin-flip and non-flip amplitudes are suppressed by ab- 
sorption similarly, because the kinematics of the data is 
rather far from the pion pole (see Table HJ. We checked 
that inclusion of absorption leaves the asymmetry well 
within the theoretical uncertainty of the calculations. 

Although the PHENIX experiment also provided data 
for An with a special selection of charged hadrons de- 
tected at pseudo-rapidities 3.0 < 77 < 3.9, we do not 
think that comparison with these data would bring new 
valuable information. Indeed, in this case summing over 
the intermediates states X one cannot rely on complete- 
ness any more, and the theoretical uncertainties would 
increase significantly. 
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VI. SUMMARY 

The results of the the PHENIX experiment at RHIC on 
single transverse spin asymmetry of neutrons produced 
in polarized pp collisions at very forward rapidities have 
brought precious information about the mechanisms of 
leading neutron production. The naturally expected pion 
pole calculated in the Born approximation does not ex- 
plain either the magnitude of the cross section, nor the 
observed azimuthal asymmetry. The absorptive correc- 
tions, which can be interpreted as the survival probabil- 
ity of a rapidity gap, substantially reduce the cross sec- 
tion improving agreement with data. At the same time, 
the relative phase shift between the spin-flip and non-flip 
amplitudes generated by absorption turns out to be far 
too small to explain the measured single transverse spin 
asymmetry. 

Another possible source of a considerable spin effect 
considered here is the interference between the ampli- 
tude of neutron production via pion and ai Reggeon ex- 
changes. Because ai has unnatural parity, it can be pro- 
duced diffractively in n+p — > ai +p, so is not suppressed 
at high cm. energy Mx- It also provides a large, close 
to maximal, relative phase shift between the non-flip oi 
and spin-flip pion exchange amplitudes. 

It turns out, however, that the oi exchange contribu- 
tion is strongly suppressed by the smallness of the diffrac- 
tive oi resonance production. Nevertheless, we found 
that it is possible to replace this resonance by np in the 
unnatural parity 1^5* state, because it forms a narrow 
peak in the Stt invariant mass distribution, so can be 
treated as an effective pole, named a, in a dispersion re- 
lation for the axial current. 

This allows to determine the a-nucleon coupling using 
PCAC, which relates the contributions of heavy states 
(saturated by the a pole and the pion pole). Additional 
information about the leptonic decay constant of a is 
obtained from the second Weinberg sum rule. 

Although the np exchange corresponds to a Regge cut. 



rather than a pole, we found its Regge intercept to be 
rather close to the one for ai Reggeon, so the phase shift 
is similar as well. 

Finally, we calculated the single transverse spin asym- 
metry at different values of the kinematic variables, s, qx 
and z, and found very good agreement with data. 
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Appendix A. FRACTIONAL NEUTRON 
MOMENTA IN PHENIX DATA 

We extracted values of z for each experimental point 
shown in Fig. [TJ from data 1], using the relation with 
the scattering angle, qt — zsm9y/l/2. The results for 
all points are presented in Table HI 

TABLE I: Values of angles of neutron detection, of neutron 
transverse momenta, and deduced values of fractional neutron 
momenta z, in the results of the PHENIX measurements [ij] 
of the neutron single transverse spin asymmetry. 



^(GeV) 


62 


200 


500 


S(mrad) 


1.89 


1.38 1.57 1.87 2.21 


1.083 1.36 1.74 2.15 


gT(GeV) 


0.05 


0.11 0.12 0.15 0.18 


0.20 0.25 0.31 0.39 


z 


0.89 


0.78 0.78 0.79 0.79 
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